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We use a combination of optical spectra, first principles calculations, and energy dependent 
magneto-optical measurements to elucidate the electronic structure and to study the phase diagram 
of NiaV^Os- We find a remarkable interplay of magnetic field and optical properties that reveals 
additional high magnetic field phases and an unexpected electronic structure which we associate 
with the strong magneto-dielectric couplings in this material over a wide energy range. Specifically, 
we observed several prominent magneto-dielectric effects that derive from changes in crystal field 
environment around Ni spine and cross-tie centers. This effect is consistent with a field-induced 
modification of local structure. Symmetry-breaking effects are also evident with temperature. We 
find M3V2O8 to be an intermediate gap, local moment band insulator. This electronic structure is 
particularly favorable for magneto-dielectric couplings, because the material is not subject to the 
spin charge separation characteristic of strongly correlated large gap Mott insulators, while at the 
same time remaining a magnetic insulator independent of the particular spin order and temperature. 
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I. INTRODUCTION 

N13V2O8 is a particularly interesting magnetic 
material i^Siiii^ both because of its unusual structure, 
which provides an example of a spin 1 system on a 
Kagome staircase lattice, and because of the rich vari- 
ety of magnetic and structural phases that are stabilized 
under different conditions. One especially interesting fea- 
ture is the occurrence of a magnetic, ferroelectric phase as 
a function of temperature and magnetic field. More gen- 
erally, coupled magnetic and electric degrees of freedom, 
flexible lattices, and magnetic frustration in multiferroics 
can result in cascades of coupled magnetic and dielectric 
transitionsi^iSiSiifliiLiSii^ The recent report^ of colossal 
low frequency (1 kHz) magneto-dielectric effects in inho- 
mogeneously mixed- valent DyM^Os is especially impor- 
tant, because it illustrates that sizable dielectric contrast 
can be achieved by physical tuning through an unusual 
commensurate-incommensurate magnetic transition and 
is facilitated by a soft lattice. The 300 K low-frequency 
magneto-dielectric effect in mixed- valent LuFe2C>4 has 
also attracted attention due to the very low magnetic 
fields needed to achieve dielectric contrast^ N13V2O8 
is another system where the temperature and field de- 
pendence of the spontaneous polarization shows a strong 
coupling between magnetic and ferroelectric order^*£ 
This coexistence is unusual and appears only when cer- 
tain symmetry conditions are fulfilled! 3 " That the effect 
can be controlled with an external magnetic field makes it 
attractive for device applications. Based upon our previ- 
ous work with inhomogeneously mixed- valent K^VsOsr^ 



the significant coupling between spin, lattice, and charge 
degrees of freedom make M3V2O8 an excellent candidate 
for discovery of higher energy magneto-dielectric effects. 

Figure QJa) shows the orthorhombic (Cmca) crystal 
structure of NiaA^Os- It consists of Kagome layers of 
edge sharing NiOg octahedra separated by nonmagnetic 
VO4 tetrahedra. N13V2O8 is considered to be a Kagome 
staircase compound due to buckling of the lattice per- 
pendicular to the a axis. There are two distinct types 
of Ni 2+ (S = 1) centers, which we refer to as "spine" 
and "cross-tie" sites. The "spine" sites run along the a 
axis. A view of the Kagome staircase showing only the Ni 
atoms is displayed in Fig. ^b). Note that the spine and 
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FIG. f: (Color online) (a) Crystal structure of Ni 3 V 2 08. (b) 
View of Kagome staircase showing only the Ni atoms. 



2 



cross-tie sites have very different focal symmetries. The 
spin ordering arrangements in M3V2O8 have been exten- 
siveiy investigated by neutron scattering and derive from 
various local and long range exchange, spin anisotropy, 
Dzyaloshinskii-Moriya, dipolar, and frustration effects. 5 
The H-T phase diagram (for H\\b) is shown in Fig. 
Hfb). Two incommensurate phases are observed beiow 
the paramagnetic phase (PM). HT1 is the (longitudinal) 
high temperature incommensurate phase, and LTI is the 
(spiral) low temperature incommensurate phase. The 
latter displays ferroelectricity. Commensurate canted an- 
tiferromagnetism is observed in the C phase. 
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FIG. 2: (Color online) (a) Specific heat of NisV^Os as a func- 
tion of temperature at H = 8 T. The arrow indicates a pro- 
posed phase transition temperature, (b) H-T phase diagram 
of Ni3V20s for H\\b with additional data points extracted 
from specific heat (solid triangles) and magneto-optical mea- 
surements (solid circles). Solid lines are discussed in the text 
and the dashed line is a guide to the eye. A new high field 
(HF) phase is also indicated. 

Here we use a combination of optical spectroscopy, first 
principles calculations, and energy dependent magneto- 
optical measurements to elucidate the electronic struc- 
ture and to study the phase diagram of M3V2O8. We 
find a remarkable interplay of magnetic field and opti- 



cal properties that reveals additional high magnetic field 
phases and an unexpected electronic structure which we 
associate with the strong magneto-dielectric couplings 
in this material over a wide energy range. Specifically, 
we observed several prominent magneto-dielectric effects 
that derive from changes in crystal field environment 
around Ni spine and cross-tie centers. This effect is con- 
sistent with a field-induced modification of local struc- 
ture. Symmetry-breaking effects are also evident with 
temperature. Even though both M3V2O8 and the pro- 
totypical Mott insulator NiO are based on Ni 2+ ions oc- 
tahcdrally coordinated by O with similar bond lengths, 
we show that the basic electronic structures of these two 
materials are very different. NiO has a large band gap 
and local moments that derive from strong Coulomb in- 
teractions and separation of spin and charge degrees of 
freedom. This separation is theoretically interesting but 
disadvantageous if one wants to promote coupling effects. 
The Slater insulator is at the opposite extreme. Here, a 
small band gap is present in the ground state, but it is 
a direct result of a specific spin ordering. In contrast, 
we find M3V2O8 to be an intermediate gap, local mo- 
ment band insulator. This electronic structure is partic- 
ularly favorable for magneto-dielectric couplings, because 
the material is not subject to the spin charge separation 
characteristic of strongly correlated large gap Mott in- 
sulators, while at the same time remaining a magnetic 
insulator independent of the particular spin order and 
temperature. 



II. METHODS 

A. Crystal Growth 

The M3V2O8 single crystal samples were prepared 
from a BaO-V20s flux. The crystals used for these mea- 
surements were grown as platelets, with their largest 
faces (a few square millimeters) perpendicular to the 
crystallographic b axis. 



B. Spectroscopic Investigations 

Near normal ac-plane reflectance of Ni3V20s was 
measured over a wide energy range (3.7 meV - 6.5 
eV) using several different spectrometers including a 
Bruker 113 V Fourier transform infrared spectrometer, a 
Bruker Equinox 55 Fourier transform infrared spectrom- 
eter equipped with an infrared microscope, and a Perkin 
Elmer Lambda 900 grating spectrometer, as described 
previously^ The spectral resolution was 2 cm" 1 in the 
far and middle-infrared and 2 nm in the near-infrared, 
visible, and near- ultraviolet. Optical conductivity was 
calculated by a Kramers-Kronig analysis of the measured 
reflectance^ An open flow cryostat provided tempera- 
ture control. 
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The magneto-optical properties of N13V2O8 were in- 
vestigated between 0.75 and 4.1 eV using a 3/4 m grat- 
ing spectrometer equipped with InGaAs and CCD detec- 
tors and a 33 T resistive magnet at the National High 
Magnetic Field Laboratory (NHMFL), in Tallahassee, 
FL. 150, 300, and 600 line/mm gratings were used, as 
appropriate. Experiments were performed with polar- 
ized light (E\\a and E\\c) in the temperature range be- 
tween 4.2 and 15 K for applied magnetic fields up to 
30 T (H || b). With light shining on the sample, we esti- 
mate that the base temperature was ~5.0 K. The field- 
induced changes in the measured reflectance were stud- 
ied by taking the ratio of the reflectance at each field 
with the reflectance at zero field, i.e., \R(H)/K(H = 
T)]. This normalized response is a sensitive way to 
view the field-induced optical changes^ Since e(E) = 
ti{E) + ie 2 (E) = e x (E) + ^-a x {E), it is clear that the 
field-induced changes in reflectance translate into finite 
energy magneto-dielectric effects. To obtain the 30 T 
optical conductivity (cti) and dielectric response (ei), we 
renormalized the zero-field absolute reflectance with the 
high-field reflectance ratios, and recalculated <j\ and e\ 
using Kramers-Kronig techniques 



C. Electronic Structure Calculations 

The calculations were done with the experimental 
crystal structural using the general potential linearized 
augmented planewave (LAPW) method, with local 
orbitalsfi2i2S as implemented in the WIEN2K program^ 
The augmented planewave plus local orbital extension 
was used for the O s and p states, the metal d states and 
the semi-core levels^ The valence states were treated in 
a scalar relativistic approximation, while the core states 
were treated relativistically. The calculations were done 
self-consistently using well converged basis sets and zone 
samplings based on 144 k-points in the irreducible wedge. 
A more dense mesh of 1200 k-points was used for the cal- 
culations of the optical conductivity. The LAPW sphere 
radii were 1.8, 1.7 and 1.5 ao for Ni, V and O, respec- 
tively. Calculations were done in the local spin density 
approximation (LSDA) and the LDA+U method. For 
the LDA+U calculations, we used values of U ranging 
from 5 eV to 7 eV on the Ni sites. The results shown are 
for U = 5 eVM. 



III. RESULTS AND DISCUSSION 

A. Optical and Electronic Properties of NisVaOs 

Figure IUa) displays the polarized optical conductivity 
of Ni 3 V 2 8 in the PM phase at 300 and 12 K. The spec- 
tra show directionally dependent vibrational and elec- 
tronic excitations with an optical energy gap of ^0.35 
eV. Based on our electronic structure calculations (dis- 
cussed in detail below) and comparison with chemically 
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FIG. 3: (Color online) (a) Polarized optical conductivity of 
N13V2O8 at 300 and 12 K, extracted from reflectance mea- 
surements (inset) by a Kramers-Kronig analysis. The energy 
range of our magneto-optical measurements is indicated by 
the arrow (0.75 - 4.1 eV). (b) A close-up view of the E\\a 
optical conductivity at 12, 8, and 5 K, corresponding to Ni 
(spine) d to d on-site excitations. Inset shows a close-up view 
of Ni (cross-tie) d to d on-site excitations at T = 300, 100, 50, 
and 12 K, respectively. These curves are offset for clarity. 



similar Ni-containing compoundsj2ai25i2£i2i we assign the 
excitations centered near ^0.75 and ^1.35 eV to Ni d 
to d on-site excitations in the minority spin channel on 
cross-tie and spine sites, respectively. These excitations 
are optically allowed due to the modest hybridization be- 
tween the Ni d and O p states. We assign the ~3.0 eV 
feature in the optical conductivity spectrum to a combi- 
nation of O 2p to Ni 3d and O 2p to V 3d charge transfer 
excitations. The ~4.4 eV feature derives from O 2p to V 
3d charge transfer excitations. 

Transition metal dtod on-site excitations are sensitive 
indicators of the local crystal field environment. Figure 
EJb) shows a close-up view of the Ni (spine) d to d on- 
site excitations in the -E||a optical conductivity at 12, 8, 
and 5 K. This structure is broad and featureless in the 
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PM phase (12 K), but it splits into at least four different 
components at 8 K (HTI phase). The splitting becomes 
more pronounced at 5 K (LTI phase). This splitting is 
due to a local structure distortion around the Ni cen- 
ter with decreasing temperature and demonstrates that 
the NiC>6 environment is different in the PM, HTI, and 
LTI phases. Although it is at the limit of our detector 
range, the trailing edge of the ~0.75 feature also dis- 
plays splitting below 12 K (not shown) Interestingly, 
the Ni (cross-tie) d to d on-site excitation also shows a 
notch-like structure that develops between 100 and 50 K 
(inset of Fig. EJb)), indicating that a weak structural dis- 
tortion in the cross-tie direction precedes the cascade of 
low temperature magnetic phases. Anomalous interme- 
diate temperature phonon shifts have also been reported 
in the RMn 2 05 (R = Bi, Eu, Dy) family of frustrated 
multiferroics ^ 



B. Electronic Structure Calculations of NisV^Os 

As mentioned above, there are two nonequivalent Ni 
sites in the unit cell, a spine site (4 atoms per cell, de- 
noted Ni2 in the following) and a cross-link site (2 atoms 
per cell, denoted Nil). As mentioned, N13V2O8 is a lo- 
cal moment magnet with a complex magnetic phase di- 
agram. To assign the peaks in the optical spectrum, we 
did calculations for two magnetic structures: (1) all spins 
on each Ni sublattice aligned ferromagnetically, but with 
the two sublattices opposite to each other (denoted FiM, 
in the following) and (2) a ferromagnetic ordering (FM). 
In the LSDA, the FiM ordering had lower energy by 7 
meV per formula unit (note that there are two formula 
units per cell). This shows an antiferromagnetic interac- 
tion between the two sublattices. 

The calculated total and projected densities of states 
(DOS) for the FiM ordering are shown in Fig. 0|and Fig. 
El within the LSDA and the LDA+U (U = 5.0 eV) ap- 
proximation, respectively. In addition to the U = 5 eV 
DOS shown, we did calculations for U = 6 eV and U = 
7 eV (see below). Insulating behavior is found in both 
approximations. The gaps are -E 9 (LSDA) = 0.30 eV and 
£ 9 (LDA + U) = 1.90 eV. There are only relatively small 
changes in the O 2p bands between these two approxima- 
tions. The O 2p valence bands extend from ~ -6.5 eV to 
the valence band edge in the LSDA, and from ~ -6.0 eV to 
the band edge in the LDA+U approximation. Wang and 
co- workers™ investigated M3V2O8 as a potential pho- 
tocatalyst for water splitting using optical spectroscopy 
from 300 nm to 850 nm, and density functional calcula- 
tions. However, they did not include spin-polarization in 
their band structure calculations, and this is needed to 
obtain the exchange splitting of the Ni d states. As a 
result, they obtained a metallic electronic structure with 
a high density of Ni d states at the Fermi energy in dis- 
agreement with experiment. 

The compound is described as Ni^Vj 8 ~ . As shown 
in the projected DOS, we find the V d bands well above 
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FIG. 4: (color online) Total (top) and projected DOS (bot- 
tom) as obtained with the LSDA for the FiM ordering (see 
text). The total DOS is on a per formula unit spin basis. The 
projections are onto the LAPW spheres, and are given per 
atom. 



valence bands in both the LSDA and LDA+U calcula- 
tions (from ~ 2 eV to ~ 5 eV, relative to the valence 
band edge), in agreement with this ionic model. In all 
calculations, we find integer spin magnetizations of 2 [ib 
per Ni ion, as expected for Ni 2+ . In the crystal field, the 
main Ni d splitting is due to the local coordination, which 
gives a lower 3-fold degenerate t 2g and a higher two-fold 
degenerate e g manifold per spin. The Ni d bands are al- 
ready very narrow in the LSDA. This is the reason for the 
clean gap. In the FM case, the band gap E g is reduced by 
only 0.11 eV relative to the FiM case. Presumably, an- 
tiferromagnetism within the individual sublattices would 
narrow the bands further and increase E g relative to the 
FiM ordering. Thus, already at the LSDA level, Ni 3 V 2 8 
is a local moment insulator. The majority spin Ni d lev- 
els are centered at -1.7 eV (t 2g ) and -0.6 eV (e g ), for both 
Ni sites. The corresponding minority spin levels are cen- 
tered at -0.4 eV (t 2g ) and 0.8 eV (e g ). 

The main difference between the LSDA and LDA+U 



electronic structure is in the position of the Ni d bands. 
The majority spin Ni manifolds are shifted deep into the 
valence band. These are centered at ~ -5.5 eV fo g ) and 
~ -2.5 eV (e g ). The minority spin states are split apart 
by the Coulomb interaction and are now centered at ~ 
-1 eV fog) and ~ 2.5 eV (e g ). For this value of U the 
conduction band edge onset derives from both V d and 
Ni d states. Calculations with U = 6 eV and U = 7 eV 
increase the differences from the LSDA, pushing the Ni d 
conduction band onset above the V d onset and driving 
the occupied Ni d manifolds deeper into the O p bands. 

The calculated LSDA and LDA+U optical spectra for 
FiM ordering are shown in Fig. EJ Besides the larger 
LDA+U band gap, the shapes of the spectra are very 
different. In the LSDA the onset has both Ni d - Ni d and 
charge transfer character. The first structure just above 1 
eV is due to the structure in the Ni d to Ni d excitations. 
It contains two peaks. The Nil (cross-link) site has a 
lower primary crystal field splitting than the Ni2 (spine) 
site, as may be seen from the projected DOS. Therefore, 
the lower energy peak comes more from the Nil and the 
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FIG. 5: (color online) Total (top) and projected DOS (bot- 
tom), as in Fig. 0] but with the LDA+U method, [7 = 5 
eV. 
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FIG. 6: Real part of the optical conductivity within the LSDA 

(top) and LDA+U, U = 5 eV (bottom), for Ni 3 V 2 8 with 

the FiM ordering. A 0.1 eV broadening was applied. Two 

definitions of the crystallographic axes are in literature. We 

use the one where a = 5.936 A, 6 = 11.420 A, and c = 8.240 
X i 



higher peak from the Ni2. The second main structure 
above 2 eV, which is weaker, is of mainly charge transfer 
character involving Ni and V, and the onset starting at 
3 eV is mainly of charge transfer character into the V 
d bands. In the LDA+U case the spectrum starts at 
~ 2 eV and is dominated throughout by charge transfer 
excitations. 

The LSDA spectrum is clearly in better accord with 
the experimental spectrum than the LDA+U spectrum, 
both in terms of the gap and in terms of the peak struc- 
ture. The main differences between LSDA and experi- 
ment are in the region of the second peak (above 2 eV) 
and in the onset of the O 2p to V d charge transfer exci- 
tations. Considering that V is in a d° configuration, the 
underestimate of this gap is not surprising and presum- 
ably reflects just the generic LDA band gap underesti- 
mation, which is seen in other d° oxides. The relatively 
good agreement of the LSDA Ni d position and the fun- 
damental gap with experiment is, however, unexpected 
and deserves comment. 

NiO, which like Ni3V20s, has Ni 2+ , six-fold coordi- 



nated by O, is a prototypical Mott insulator. LSDA cal- 
culations for antifcrromagnetic NiO show a small band 
gap and features in the valence band dispersions in ac- 
cord with experimental However, the gap obtained is 
much smaller than the experimental gap of ~ 4 eV, 
and furthermore it is dependent on the specific mag- 
netic ordering. Furthermore, the LSDA gap of NiO is 
of incorrect d-d character rather than the charge trans- 
fer character observed in experiment. This shows an 
essential role for beyond LSDA Coulomb correlations 
in describing the electronic structure of NiO. Various 
methods, including the LDA+U approximation, for in- 
corporating these have been developed and tested us- 
ing NiOi^VV^r^i^VVi As mentioned, NiO and 
M3V2O8 both have Ni 2+ in approximately octahedral 
O cages. In NiO, the Ni band width is larger than in 
M3V2O8, which explains why the LSDA is able to pro- 
duce local moment insulating behavior in this compound 
but not in NiO. The narrower bands in Ni 3 V20s are re- 
lated to the bonding topology. In NiO, the O are six fold 
coordinated with 90° and straight 180° bonds. Ni 3 V 2 8 
has lower O coordination and bent bonds. In metals, the 
criterion for dividing strongly correlated materials from 
band metals involves U 8 /W, where W is the band width 
and U s is the screened Coulomb potential, where the "s" 
is simply to distinguish this parameter from the U used 
in the LDA+U approximation. In general, Mott insula- 
tors derived from metals with U s /W cannot be described 
in the LSDA and require the use of the LDA+U approx- 
imation or more sophisticated approaches, while materi- 
als with weaker correlations are often better described in 
LSDA than in LDA+U calculation. Since the band width 
of NiaV208 is not larger than that of NiO, the implica- 
tion of our results is that the effective screened Coulomb 
interaction U s on at least one of the Ni sites is sufficiently 
lower than in NiO to cross-over from strongly correlated 
Mott insulator to band insulator. We speculate that this 
could be due to the lower cation coordination of 2 ~ in 
M3V2O8. Specifically, in NiO, six Ni atoms compete for 
bonding with each O (two per O p orbital). The average 
Ni-0 bond length in M3V2O8 is ~ 0.03 A shorter than 
in NiO, consistent with this from the bond valence point 
of view 42 This is not the case in M3V2O8, which may 
enable the O cage to more effectively screen the bare U. 
As mentioned, we find M3V2O8 to be a local moment 
insulator with weak sensitivity of the band structure to 
the magnetic order at the LDA level. 

Figure [7| shows the calculated optical spectra for FM 
ordering. As may be seen it is qualitatively similar in 
structure to that of FiM ordering, but shows noticeable 
quantitative differences. Most notably there are changes 
in the shape of the low energy Ni d - d peak, a suppres- 
sion of the second peak and changes in the shape of the 
onset of the higher energy charge transfer onset. Since 
neither the FiM or the FM ordering are the ground state 
ordering, these differences should be viewed as indicating 
the errors on our calculation when comparing with the 
experimental zero field spectrum. They also support the 
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FIG. 7: Real part of the optical conductivity within the 
LSDA, as in Fig. 03 but with FM ordering. 

experimental finding of substantial magnetochromic ef- 
fects in this material. Interestingly, we find considerable 
orientation dependence of the magnetochromic effects in 
this material. 

Optical excitations are charge excitations. Since Hub- 
bard correlations generally harden charge degrees of free- 
dom leading to spin charge separation one may expect 
stronger charge spin coupling, and perhaps stronger mag- 
netochromic effects in less correlated materials. This is 
also relevant to the magnetoelectric coupling observed in 
Ni 3 V 2 08, since strong Coulomb correlations, if present, 
would reduce the spin-lattice coupling, so that while fer- 
roelectricity would still be induced by the incommensu- 
rate magnetic order, it would be weaker due to reduced 
coupling. Based on our results, we find that the Ni 
sites in Ni3V2 0s are not subject to strong correlations, 
and this underlies the substantial magnetochromic effects 
seen in this material. 



C. Magneto-Optical Properties of NI3V2O8 

Figure [Sfa) and (b) shows the reflectance ratio of 
Ni 3 V 2 8 , R(H)/R(H=0 T), at 5 K. Rich field-induced 
changes in reflectance are observed over the full energy 
range in both polarizations. The most significant modifi- 
cations are near 0.75 and 1.35 eV, providing a preview of 
the physical origin of the field-induced effects. Based on 
the aforementioned positions of the magneto-optical fea- 
tures, we attribute these changes to field-induced modi- 
fications of Ni (cross-tie) d to d and Ni (spine) d to d on- 
site excitations. At 30 T, the reflectance deviates from 
unity by ±15% near 1.35 eV. On the other hand, re- 
flectance changes in the visible range (centered at ~2.2, 
2.6, and 3.3 eV) are much broader and more modest in 
size (up to ~4% for 30 T at 5 K), similar in magnitude 
to "magnetochromic effects" reported on other complex 
materials ^SSiiiii These features can be attributed to 
field-induced changes in O 2p to Ni 3d and O 2p to V 3d 
charge transfer excitations. 
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FIG. 8: (Color online) The normalized magneto-optical re- 
sponse, K(H) /K(H—0 T), of NiaV208 in an applied mag- 
netic field (H \\b) from to 30 T at 5 K (a) for light polarized 
along the a direction, (b) for light polarized along the c direc- 
tion. The insets show close-up views of the magneto-optical 
response near the Ni (cross-tie) and Ni (spine) d to d on-site 
excitations. Data were taken in 2 T steps, but only represen- 
tative curves are displayed for clarity. 



It is desirable to correlate field-induced reflectance ra- 
tio changes with the dispersive and lossy response of the 
material. Figure |§1 displays the polarized optical con- 
ductivity of Ni 3 V 2 8 at 5 K at and 30 T (H\\b), ex- 
tracted from a combination of absolute reflectance mea- 
surements, the reflectance ratio data of Fig. [SJa) and (b), 
and a Kramers-Kronig analysis and calculation of the op- 
tical constants. From the optical conductivity, we can 
immediately confirm that the aforementioned magneto- 
optical effects correspond to field-induced modifications 
of the Ni (cross-tie) d to d, Ni (spine) d to d on-site ex- 
citations, and O 2p to Ni 3d and O 2p to V 3d charge 
transfer excitations. The insets of Fig. [5] show close- 
up views of the ~1.35 eV Ni on-site excitation splittings 
that point toward a field-induced local distortion of the 
NiC>6 octahedral As mentioned previously, these Ni d 
to d excitations are allowed because oxygen hybridization 
modifies the matrix elements. The charge transfer exci- 
tation centered at ~3.0 eV also broadens and blueshifts 
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FIG. 9: (Color online) Polarized optical conductivity of 
Ni 3 V 2 8 at 5 K for H = and 30 T (H\\b). E\\a and E\\c 
polarizations are indicated by dashed and solid lines, respec- 
tively. The upper inset shows a close-up view of the energy 
region near the Ni (spine) d to d on-site excitation for H = 
and 30 T. The lower inset shows a close-up view of the Ni 
(spine) d to d on-site excitation for T = 12 and 5 K at 30 T. 



with applied field, and a small "notch" develops at ~ 2.2 
eV. The reflectance ratio trends discussed above also cap- 
ture important dispersive effects. Figure lTUT a) shows the 
high energy magneto-dielectric response of N13V2O8 at 5 
K for E\\a. Although dielectric contrast (Aei/ei) is sub- 
tle in the majority of the visible region (Fig. llUf c)). the 
applied field changes the dielectric constant significantly 
near the Ni d to d on-site excitations (inset, Fig. HOf a). 
and Fig. HUf b)). At some energies, dielectric contrast is 
on the order of 15% (Fig. HUf b)). These results demon- 
strate that magneto-dielectric effects are not confined to 
kHz frequencies and that the response may be tunable. 
Both aspects may be useful for device applications (for 
instance, magnetically controlled ferroelectric memory). 
This remarkable interplay between the magnetic field and 
optical constants of M3V2O8 seems to be facilitated by 
the fact that this material is in the intermediate coupling 
regime, as discussed below. 

Comparison of the reflectance ratio data in the upper 
insets of Fig. |5fa) and (b) provides insight into possi- 
ble magnetoeleastic coupling in N13V2O8. The structure 
centered at ~0.83 eV is large in the c (cross-tie) direc- 
tion, whereas it is much smaller along the a (spine) direc- 
tion. It grows quickly with applied field. The splitting of 
this structure (~12 meV (97 cm -1 )) corresponds to a vi- 
brational energy scale. Candidates for possible coupling 
might therefore include the displacement of Ni (cross- 
tie) atoms or a low-frequency O-Ni-0 bending mode. In 
the presence of a soft lattice, magnetoelastic coupling 
may lead to a local distortion of the NiOe octahedra. 
Note that the PM phase of a frustrated magnetic sys- 
tem often displays significant spin correlations^ The 
PM phase in M3V2O8 may be similar, consistent with 
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FIG. 10: (Color online) (a) Dielectric constant of M3V2O8 
at 5 K for light polarized along the a direction for H = 
and 30 T (H\\b). The inset shows a close-up view of the 
energy region near the Ni (spine) d to d on-site excitation. 
Close-up views of the high-energy dielectric contrast, Aei/ei 
= [ei(E,H) - ei(E,0)]/ei(E,0), near (b) the Ni (spine) d to 
d on-site excitation and (c) the charge transfer excitations. 



the presence of magnetoelastic coupling. Direct mea- 
surements of high field vibrational properties of NI3V2O8 
are clearly desirable. We note that spin-lattice coupling 
leads to structural distortions in the frustrated magnetic 



system ZnCr 2 0, 



45.46 



Magnetoelastic coupling also plays 



an important role in geometrically frustrated magnetic 
systems such as HoMnOa and RM^Os (R = Tb, Ho, 

p..j 7.8.10.29.47.48 

Spin-charge coupling in N13V2O8 can be further in- 
vestigated by quantifying the subtle changes in the Ni 
(spine) d to d on-site excitation with applied field. Fig- 
urc llll shows a close-up view of polarized optical conduc- 
tivity of Ni 3 V 2 8 at 5 K for H = 0, 10, 20, and 30 T 
(H\\b). The effect of applied field is non-monotonic and 
correlates with changes in magnetic order. Peak fits, us- 
ing four model oscillators, were used to elucidate these 
trends (inset, Fig. Ill|) . The dependence of peak po- 
sitions on the applied field clearly demonstrates an in- 
terplay between electronic and magnetic properties. For 
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FIG. 11: (Color online) Polarized optical conductivity of 
N13V2O8 at 5 K for H = 0, 10, 20, and 30 T (H\\b). The 
inset shows the peak position of fitted oscillators as a func- 
tion of applied magnetic field. The shaded regions represent 
transition fields, from which the extended H-T phase diagram 
in Fig. HJb) was generated. 



instance, at 5 K, the data shows discontinuities near 7 
and 15 T, indicated by the grey shaded regions in the 
inset. These transition regions move to higher fields with 
increasing temperature. We note that the low field phase 
boundaries determined in this manner are in good overall 
agreement with the previously reported thermodynamic 
data. This optical properties work also allows us to fol- 
low the C/LTI phase boundary to higher fields, where it 
crosses over both the HTI and PM phases. Finally, these 
studies also identify a new "High Field" (HF) region, 
nestled between C and the cascade of LTI, HTI, and PM 
phases at higher temperature (Fig. I^b)). Reexamina- 
tion of previous specific heat measurements on N13V2O8 
are consistent with the presence of an additional weak 
transition for H = 8 T, as shown by the arrow in (Fig. 
IUa)). A similar transition was also observed at H = 7 
T. 

As discussed in the Introduction, Ni 3 V 2 08 has a com- 
plex H-T phase diagram with several low field phases 
arising from the complex spin structures!^ With the new 
transition fields extracted from the magneto-optical data 
at 5, 8, and 12 K (solid circles) and specific heat data 
(solid triangles) , we extended the H-T phase diagram of 
Ref. 4] for H\\b, as shown in Fig. H£b). We used the 
combined optical and specific heat data to establish the 
boundaries for the newly identified high field phase. It 
has been predicted^ that the phase boundary for the LTI 
phase should be quadratic in H, so that (T c (H)-T c (0)) 
~ H 2 . Rather than this quadratic dependence, we find 
that this phase boundary can be fit by (T c (H)-T c (0)) ~ 
H 2 - 5 . This fit is shown by the lower solid line in Fig. 
|2Jb). The origins for this discrepancy are unclear at this 
time. The upper solid line in Fig. |2{b) is a guide to the 



eye for the high field (HF) phase boundary. To the best 
of our knowledge, there has been no discussion of any 
possible functional form that might apply to this bound- 
ary. Combined magnetization and specific heat measure- 
ments in high magnetic fields will be required to precisely 
determine the magnetic field dependence of both phase 
boundary lines, which in turn will provide mechanistic 
information on the underlying magnetic transitions. Ad- 
dition of the new data points from optical and specific 
heat studies also demonstrates the existence of a new 
HF phase that overlays the other low field regions. Neu- 
tron experiments are needed to assess magnetic order in 
the HF phase. 

In addition to the 5 K (LTI phase) data discussed in 
detail here, we also collected full data sets in the HTI 
and PM phases of Nis^V^Os- Strong magneto-optical ef- 
fects were observed at both 8 and 12 K. As shown in 
Fig. the magneto-optical effects in the PM phase 
are qualitatively similar to those of the LTI and HTI 
(not shown) phases, although the overall size of the field- 
induced modification is smaller. This result demonstrates 
that an applied field induces a local structural distortion 
around the Ni sites in all of the low temperature phases. 
In otherwords, the field-induced distortion is not unique 
to the LTI (ferroelectric) phase with spontaneous polar- 
ization. To support this claim, we refer the reader to 
the inset of Fig. 112( b)). where the zero-field ^1.35 eV 
excitation is featureless. An applied field splits this struc- 
ture into several distinct components, demonstrating the 
field-induced distortion of the Ni (spine) crystal field en- 
vironment even in the PM phase. 



IV. CONCLUSION 

We investigated the optical properties, electronic 
structure, and energy dependent magneto-optical re- 
sponse of Ni3V2C>8 to elucidate the electronic structure 
and to study the phase diagram. The spectra exhibit 
features centered at ~0.75 and ~1.35 eV that we assign 
as Ni (cross-tie and spine) d to d on-site excitations in 
the minority spin channel. O 2p to Ni 3d and O 2p to 
V 3d charge transfer excitations appear at higher en- 
ergy. Extensive analysis of splitting patterns of the Ni 
(spine) d to d excitation in the PM, HTI, and LTI phases 
demonstrates that the local Ni environment is sensitive 
to magnetic order even at zero field. A splitting of Ni 
(cross-tie) d to d on-site excitation between 50 and 100 
K also points toward a weak local structural distortion 
that precedes the cascade of low temperature magnetic 
phases. Although both N13V2O8 and the prototypical 
Mott insulator NiO are based on Ni 2+ ions octahedrally 
coordinated by O with similar bond lengths, the basic 
electronic structures of these two materials are very dif- 
ferent. NiO has a large band gap and local moments that 
derive from strong Coulomb interactions and separation 
of spin and charge degrees of freedom, while, in contrast, 
we found M3V2O8 to be an intermediate gap, local mo- 
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FIG. 12: (Color online) (a) The normalized magneto-optical 
response, R(H)/~R(H=0 T), of M3V2O8 in an applied mag- 
netic field (H\\b) from to 30 T at 12 K for light polarized 
along the a direction, (b) Polarized optical conductivity of 
Ni 3 V 2 8 at 12 K for H = and 30 T (H\\b). Inset shows 
a close-up view of the magneto-optical response near the Ni 
(spine) d to d on-site excitation. 



ment band insulator. This electronic structure is partic- 
ularly favorable for magneto-dielectric coupling, because 
the material is not subject to the spin charge separation 
of the strongly correlated large gap Mott insulator, while 
at the same time, remaining a magnetic insulator inde- 
pendent of the particular spin order and temperature. 
The remarkable interplay of magnetic field and optical 
properties in Ni3V20s reveals an additional high field 
phase and an unexpected electronic structure, which we 
associate with the strong magnetodielectric couplings in 
this material. We discovered several prominent magneto- 
optical effects that derive from changes in crystal field 
environment around Ni (spine and cross-tie) centers due 
to a field- induced modification of local structure of NiOe- 
The magnetoelastic mechanism, responsible for the field- 
induced distortion of the NiOg building block unit, is 
active in the paramagnetic phase as well. 
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